Quasi-Thermal Neutrinos from Rotating Proto-Neutron Stars 
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Rotating and magnetized proto-neutron stars (PNSs) may have neutron-loaded relativistic winds. 
These neutrons should cause inelastic collisions in the stellar material or around the termination 
shock, producing ~ 0.1 — 1 GeV neutrinos, without relying on cosmic-ray acceleration mechanisms. 
We show that PINGU and Hyper-Kamiokande can detect such neutrinos from nearby supernovae, 
by reducing the atmospheric neutrino background via coincident detection of MeV neutrinos or 
gravitational waves and optical follow-up observations. Detection of these quasi-thermal neutrinos 
would provide important clues to the physics of magnetic acceleration, nucleosynthesis, the relation 
between supernovae and gamma-ray bursts, and the properties of newly born neutron stars. 
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Proto-neutron stars (PNSs) are produced in core- 
collapse supernovae (SNe), and cool via radiation of 
MeV neutrinos on a timescale of ~ 10 — 100 s [e.g., 
G], A fraction of these thermal neutrinos deposit 
their energy in the PNS atmosphere, driving a wind with 
mass loss rate M that injects energy into the shocked 
stellar material and forms a wind-driven bubble in the 
SN cavity [e.g.,[3|- For non-rotating and non-magnetic 
PNSs, the wind kinetic energy is tiny compared to the 
SN explosion energy, and the wind is non-relativistic 
throughout the cooling epoch [e.g., |4|. However, if 
PNSs are rotating and magnetized, they transition from 
non-relativistic and thermally-driven to relativistic and 
Poynting-dominated [e.g., |5j, |6j. The magnetization of 
the flow is @, H [H 
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where R ns is the PNS radius, £?dip is the surface dipolar 
field strength, f2 is the angular frequency of the PNS, 
and / op takes into account that the outflow comes from 
only the open fraction of the PNS surface. The transition 
from non-relativistic to relativistic occurs at a ~ 1 where 
the Alfven speed becomes ~casM decreases. 

Rapidly rotating and/or strongly magnetized PNSs 
may be related to various explosion phenomena. Esti- 
mates suggest that > 10% of SNe lead to magnetars with 
£>di P ~ 10 14 ~ 15 G [7j, which may be generated by rapid 
rotation via the dynamo mechanism [8J. Rotation and/or 
magnetization may also modify the explosion dynamics. 
If the wind power exceeds ~ 10 48 erg s _1 , the wind can be 
collimatcd , and these jets may lead to gamma-ray 

bursts (GRBs)[l, H [ll[ . If the rotation rate or magnetic 
field strength is not sufficiently high, a quasi-spherical 
wind or a pair of choked jets hidden by the stellar ma- 
terial may result, potentially leading to hypernovae and 
super-luminous SNe [1, @, [l3| . 

Neutrons are loaded in the PNS wind by neutrino heat- 



ing, and they may also be important for synthesis of nu- 
clei [14( . We show that relativistic neutrons are produced 
by magnetic acceleration of the flow as a becomes ^> 1, 
and that ~ 0.1 — 1 GeV neutrinos should be generated 
via the np reaction. These neutrinos can be detected 
by planned facilities such as PINGU [HI and Hyper- 
Kamiokande (HK) [13], and their characterization would 
allow us to probe the physics of PNSs, the early wind/jet 
dynamics, and magnetic acceleration in an environment 
inaccessible to photons. Throughout this work, we use 
Q x = Q/10 x in CGS units unless otherwise specified. 

Relativistic PNS winds. — Mass-loss occurs during 
the PNS cooling phase by neutrino heating, mainly via 
v e n f± e~p and v e p ^± e + n. In unmagnetized winds @,@| 
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where L v is the neutrino {v e + v e ) luminosity and e v is the 
typical neutrino energy, which can be affected by rotation 
and magnetization [17j. Additional heating due to in- 
elastic electron scattering gives a correction / cs = 1 + e cs . 
The mass-loss rate is reduced by . f nn . while enhanced by 
/ cen due to the centrifugal force [18| . We describe the 
baryon mass-loss rate by M& w M = M„/ op / cen [f|. 
By setting cr ~ 1 in Eq. (1), we find that the tran- 
sition to relativistic flow occurs when Mb w M tr ~ 



7.4 x KT b M s" 1 B. 
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The neutrino luminosity decreases gradually as a power 
law, until the PNS becomes transparent to neutrinos [3, 
6|. Once i t hin ~ 10—100 s, L v , e v and M v should decline 
rapidly. Eventually, M cannot be described by neutrino 
heating but by the Goldreich- Julian density Mqj ~ 2.5 x 
1O~ 17 M s- 1 M±,6-Bdip,i 5 Pr 2 2 i^ S! 6 0, where /x± is the 
pair-multiplicity. Note that M& = AlQjm p /(fi±m e Y e ), 
where Y e is the electron fraction. 

For the purpose of simple estimates, we hereafter as- 
sume M ns = lAM m , i?ns = 10 6 cm, L v oc r~ 1 e -i /* tw » 
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£„ oc L^ 4 @, Then, the transition time tt x , when the 
PNS wind becomes relativistic (<r ~ 1) is 



~ 3 2 S B" 4/5 P 4 / 5 f"2/5f2/5 
O.z b ^>dip, 15^-2 ./op Jccn 

x L 2 / 3 52 ( £ „ /15MeV) 4/3 / c 2 / 3 , 
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where L„o and e^o are defined at 1 s. Therefore, if mag- 
netic fields are strong and/or rotation is rapid, the PNS 
wind becomes relativistic at t <C tthin. If not, it will 
become so at ~ £ t hm as M rapidly declines. 
The a parameter at tthin is 



a(t thin ) ~ iioos^ iPil5 pr 2 7op,-i/< 

x (e u0 /l5 MeV)- 10/3 / ffi 5/3 , 
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and M(t mn ) ~ 6.5 x KT 11 M Q s" 1 f op ,-iUnL%% 

(f^o/15 MeV) 10/3 /c/ 3 for tthin = 50 s. We focus on the 
epoch t < tthin, since the flux of relativistic neutrons 
producing quasi-thermal neutrinos decreases strongly for 
t S> tthin even though the wind becomes more relativistic. 

Wind acceleration. — At t > t tI , the PNS wind is 
magnetically accelerated. It is known that magnetic en- 
ergy is not efficiently converted to kinetic energy in ideal 
MHD [20] . However, efficient acceleration is required to 
explain jets of GRBs and active galaxies, and may be 
achieved by time-dependence or magnetic dissipation [see 
a review[2l|, and references therein] . When the maximum 
Lorentz factor r max « a is achieved at the saturation ra- 
dius i? sat ~ -Rmag, we can parametrize T(r) as L 22j 



r(r) w cr 
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5.0 x 
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ter characterizing magnetic reconnection |22|, [23J . Note 
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that more efficient dissipation or radiative acceleration 
leads to indices larger than 1/3 and smaller i? sa t- Hence, 
our results below are relatively conservative. 

In our case, the wind/jet-driven bubble terminates 
at the reverse shock caused by the interaction between 
the flow and shocked/pre-shocked stellar material, which 
may be nearly-spherical or biconical 0, HI| ■ When the 
wind radius i?„ 



is smaller than i? m ag, the final T is 



T(R W ) 



13rT l / 3 P~ i/l V / ' J R 1 

100 3 r -% t rcc,-2 rl ji',10- 



-1/3 1/3 



?l/3 



(6) 



The maximum Lorentz factor is achieved only if there 
is no boundary to the flow (e.g., high wind power and 
successful jets) 0,(13] or acceleration is much faster (e.g., 
choked jets in the fireball scenario for GRBs). 

Neutrons and ions have the same outflow velocity as 
long as they are coupled with < a c \v >~ a np c [2^]. Here 
<7 c i is the elastic cross section and a np ~ 3 x 10~ 26 cm 2 is 
the inelastic cross section. However, if the neutrons are 



decoupled during the flow acceleration, then they will 
have smaller Lorentz factor than the ions. Using the 
nucleon density n w = M b / '{Airr 2 cm p ) and Eq. (5), r np w 
0.5n u) tT np (i? c iec/r) = 1 gives the decoupling radius 

R dcc ~ 6.8 x 10 s cm a^ 1 B^PZl'' f%U^JU (?) 
For i?dec < -Rujj-Rmag, the neutron Lorentz factor F„ is 



-3/4,1/2 1/4 



(8) 



Because R mag oc a 2 , it should increase rapidly in time 
(Eqs. 1 and 2). R w also increases, but more slowly. 
Hence, if P ma g(ttr) < R w (ttr) at t tr , then i? mag even- 
tually overtakes R w . On the other hand, i?dec decreases 
with time, so it crosses i? mag and R w as long as i?dec(ttr) 
is large enough (that is satisfied in our cases). We define 
t a by Rdcc{t a ) = Rmngita) and t b by Rdcdh) = Rw{h), 
respectively (see Fig. 1). (Note that t a < i t hm and 
tb ^ ^thin since a abruptly increases around tthin-) 

Neutrons and ions achieve the same final Lorentz factor 
in the early phase (t < max \t a , tb]), whereas neutrons 
have lower final Lorentz factor in the later phase because 
of decoupling. If t a < tb, we have r„ ps T(R w ) at t < 
tb and r n w r(i?dcc) at t b < t. If t b < t a , we obtain 

r„ « r(i? mag ) &tt<t a and r n « r(i? dcc ) at t a < t. if 

t a < tthin, we have 



16 s B, 



-3/5 pl/2,-1/5 
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using Eqs. (1) and (2). If t& <C tthin, assuming the typical 
velocity V « R w /t, tb is estimated to be 

U ~ 40 s B~ 1/T P 3 / 14 /- 1 / 7 f2/7,-i/i4 

lb ~ W » -D d ip, 15-2 J op, -1 Jccn e rec,— 2 

x CoS^o/lSMeV) 20 / 21 /^ 21 ^ 7 . (10) 

Neutrons decay with proper lifetime of 886.7 s, so the 
decay radius Rp « 2.7 x 10 14 cm r„.i is much longer 
than i?dcc, -Rmag and i?™, in which we are interested. 

Neutrons from disintegrated nuclei. — Neutrons 
may be produced around R w via disintegration of nuclei, 
as well as those not used up in nucleosynthesis. Fe-group 
elements are mainly created if Y e o > 0.5, while the neu- 
tron capture channel is dominant and A > 56 nuclei are 
synthesized if Y e o < 0.5. They could be seeds of ultra- 
high-energy cosmic rays (UHECRs) if they could survive 
(e.g., after successful jets break out) 

Nuclei may be spalled, which is relevant at R w < -Rdcc 
since shocked nuclei collide with neutrons. They can also 
be disintegrated by photons mainly via the giant-dipole 
resonance A-y -> A'N if 2re > 8.5 MeV (^/56) _1/6 , 
where e is the target photon energy in the stellar 
frame 26]. The wind/ jet- driven bubble may be filled 
with photons provided by the SN explosion, since 
the SN shock or radioactive nuclei provide x and/or 
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t a < t b (spherical wind/choked jets) t a > t b (successful jets/choked jets) 
R t : R 




decoupling 




FIG. 1: Schematic picture of -R mag oc a 2 , Rdcc oc a 1 and 
R w oc t 1+x . 

gamma rays, cascading down to optical photons as a 
result of thermalization. The SN temperature of w 
11 kcV £ph 4 49^ 3 /o 4 allows nuclei to survive. However, 
disintegration can be caused by high-energy photons that 
are generated via shock dissipation or magnetic recon- 
nection. We here consider effects of high-energy elec- 
trons produced around R w . The PNS wind is Poynting- 
dominated after ftr, and the magnetic field is B ~ 2.1 x 
10 8 G Sdip,i5-f > r 2 1 /o P ,-i-R^ 1 io f rom the shock-jump condi- 
tions. The reverse shock is hydrodynamically weak since 
the thermal energy in the downstream is much smaller 
than the magnetic energy, but small dissipation can be 
enough. The random energy per particle is 0.75Tm p c 2 
when ions are mainly protons, and the injection Lorentz 
factor of electrons is jej ~ 0.75e e T(m p /m e )Y e ^ 1 . 
The characteristic synchrotron energy is £ syn ~ 

4.6 MeV ee > -i^ 2 ri-B d ip,i 5 Pr 2 1 /op -i-R^^io £ ™eC 2 , so 
synchrotron cascades will be developed and nuclei typi- 
cally interact with boosted ~ m e c 2 photons. Then, using 
n A a Al « 1.4 x 1CT 27 cm 2 (A/56) 1/6 [26], the photodis- 
intcgration efficiency, /a 7 ~ hag Ay(Fn^)(R w /T 3 ), is 

f Al ~ 3(^/56) 1 /^ e rrV 3 - 1 sg )1B pr 2 3/ V^-ii^,\^i) 

Note that the Thomson opacity is tt ~ n e aT(R w /^) — 
OMYeT^R-^iMb/lO- 10 - 5 M Q s" 1 ) < a few, so par- 
ticles may be accelerated and high-energy photons can 
reach the upstream before they are thermalized. 

Eq. (11) suggests that synchrotron photons can 
disintegrate nuclei and supply neutrons with T « 
min[r(i?„,), r(i? mag )], giving an additional contribution 
of neutrinos. Hence, the fraction of neutrons (X n ) de- 
pends on outflow dynamics as well as the initial entropy 
and the initial electron fraction Y e o. Note that Jaj de- 
clines with time, and nuclei would survive at t 3> ithin- 

Neutrino production. — While the wind/jet exca- 
vates the stellar material, ions are quickly decelerated 



at the shock via radiation or collisionless processes. On 
the other hand, relativistic neutrons should be deceler- 
ated via np collisions, leading to neutrino production. 



Note that their deceleration scale 



l/(n w Q.5a r . 



longer than ~ \/{n e OT) and any relevant plasma scales. 
Even above i?dec, the neutron beam goes through the 
shocked wind, but it is damped by the stellar mate- 
rial (with mass M sn ). Neutrons leaving the flow are 
depleted until the SN ejecta reaches the pionosphere at 
« 6.5 x 10 15 cm (Msn/IOM©) 172 . 

Half of the neutron kinetic energy goes to neutrinos, 
so quasi-thermal neutrinos have a luminosity of 



Cu ~ 0.5(r„ - l)X n M b c 2 



(12) 



For t a < tf,, assuming V oc t x with small \, we have 
C v oc £-( 4 -X>/3 until t b , so C v and £ v typically decreases 
as time. Inelastic collisions occur only when r„ > 1.37, 
almost after the PNS wind becomes relativistic at ~ t tI . 
(Of course, both R w and i? mag have to be larger than 
R\ c = c/fi.) Hence, the total energy of neutrinos is £ v ~ 
0.5X„(0.37 Mb)\t tr c 2 t tI , and we have 

C° ~ 7.8 x 10 48 erg X n B«£ \ 15 PlT f%* ft£ 

x L%% 2 ( Ev0 /15 MeV^f^f- 1 , (13) 

where fb is the beaming factor and the typical neutrino 
energy is E v w 0.25to t c 2 ~ 35 MeV. If the PNS wind 
is collimated as in GRB jets, the observed fluence is en- 
hanced by f b ~ ~ 100, where such neutrinos could be seen 
as a tail by water- Cherenkov neutrino telescopes such as 
Super-Kamiokande and a next-generation detector HK. 

Since T n increases with time, higher-energy neutrinos 
are produced later, and these are more easily detected. 
However, L v and C v decline abruptly after tthin- For neu- 
trinos originating from intrinsic neutrons, the highest- 
energy contribution is generated at ~ tb, at which we 
have £ v ~ 0.5X n (T n M b )\t b c 2 tb and Eq. (10) leads to 

£™ ~ 2.5 x 10 47 erg X n B^ 15 PI^ U C_J^U 

x £Sg(Wl5 MeV) 40 / 21 ^ / 21 ^/ 7 /,- 1 . (14) 

Though Eq. (14) is valid for neutrons produced via pho- 
todisintegration of nuclei, they achieve higher Lorentz 
factors at t > % and their contribution at ~ ithin is still 
significant. Using £ v ~ 0.5X n (r(R w )M b )\t thin c 2 t thin , 



C° ~ 3.8 x 10 46 erg X n B*£^2 fiLUf^^l-* 
x < 5 9 2 (W15 McVf ^ f^R^f^, (15) 

is obtained for this case (where Uhin = 50 s is used). 
For t b < t a , C v oc r„Mfc oc t until t a and T n satu- 
= i?doc- Then, for both the origins of 



rates when R„ 



neutrons, we have £ v ~ 0.5X n (T n Mb)\t a c 2 t a and 

£ iso ~ 1 x 10 48 ersrX B 7/5 P~ 3/2 f 9/5 f 2 / 5 f 1/10 
c v ~ i.u x iu erg ^„i3 dip 15 _r_ 2 / opj _x J ccn t r0 c,-2 

x L% 3 52 (e„ /15 MeVf'fTfb 1 - (16) 
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In either of Eqs. (14)-(16), the neutrino spectrum may 
extend to ~ r n> i GeV with the typical energy E v m 
0.05r n m n c 2 ~ 6.47 GeV T„ 4 for F n » 1 0. Such 
GeV neutrinos are good targets for PINGU, a planned 
low-energy extension of IceCube, as well as HK. 

In the standard pulsar phase (t ^> tthin), M « Mgj 
and we do not expect intrinsic neutrons. The wind will 
be dominated by electron-positron pairs, though nuclei 
may be stripped from the surface. Neutrons may be sup- 
plied by photodisintegration. However, since the ambient 
photon density is already low, photodisintegration and 
neutrino production are likely to be inefficient. 

Neutrino detection. — From Eqs. (13)-(16), the 
energy fluence of quasi-thermal neutrinos per flavor is 
E 2 ^ ~ 28 erg cm" 2 kpc)- 2 . 

PINGU has sensitivity to 1 — 10 GeV neutrinos with 
an effective area of ~ 4 x 10~ 3 cm 2 for v e + v e and ~ 
2 x 1CP 3 cm 2 for + P^, respectively, at ~ 1 GeV. So, 
~ 100 £„%$ events are expected for a SN at 10 kpc. In the 
case of choked jets the energy fluence may be enhanced by 
~ 100, and the detection of GeV neutrinos becomes 
possible even for extragalactic SNe out to ~ 1 Mpc. 

HK has the fiducial volume of 0.56 Mt, so the effective 
numbers of free protons and bounded nucleons in oxygens 
are 3.7 x 10 34 and 3.0 x 10 35 , respectively. The neutrino- 
nucleon cross section for the charged-current interaction 
at 1 GeV is ~ 0.6 x 10 -38 cm 2 (averaged over v and P), 
so the effective area is ~ 2 x 10~ 3 cm 2 . Hence, we may 
expect ~ 70 £ ™\ % events for a SN at 10 kpc. In addition, 
HK could also allow us to see ~ 10 — 100 MeV neutrinos, 
through the P e p — > e + n channel. But the detection would 
be more difficult due to smaller cross sections at lower 
energies and the signal may be buried in the exponential 
tail of thermal MeV neutrinos from the PNS. 

To see the signal, it is crucial to reduce backgrounds us- 
ing space/time-coincidence. The obvious background is 
the atmospheric neutrino background (ANB). The ANB 



at GeV is « 1.3 x 10^ GeV cm" 2 s" 



for v„ 



and w 2.6 x 10 -2 GeV cm~ 2 s^ 1 sr _1 for + P^, 
respectively [23] ■ We may take the time window of 
hhin ~ 10—100 s after the explosion time that is measur- 
able with MeV neutrinos and/or gravitational waves [28j |. 
The localization is possible by optical follow-up observa- 
tions, and the angular resolution of PINGU is expected 
to be - 2 - 20 deg [l]|. Then, the ANB in this time and 
angular window, ~ 2 x 10 -3 erg cm -2 , is negligible. 

Implications and discussions. — Neutron-loaded 
relativistic winds emanate from rotating and magnetized 
PNSs during their ^ 10 — 100 s cooling epoch. Rela- 
tivistic neutrons produce ~ 0.1 — 1 GeV neutrinos via 
the np reaction as the wind/jet-driven bubble interacts 
with the surrounding stellar material. For PNSs with 
-Sdip ^ 10 14 G and P ~ 1 — 10 ms, future neutrino tele- 
scopes such as PINGU and HK may detect ~ 10 — 100 
neutrino events of energy ~ 0.1 — 1 GeV from the next 
Galactic SN. Even extragalactic SNe may be detected 



if the wind is collimated, forming choked jets. Al- 
though not all SNe are born with such high B^ip and 
short P, > 10% of SNe may be [7], and we may ex- 
pect 0.02 — 0.05 yr _1 for the birth of magnetars within 
5 Mpc [29] . In addition, given the ANB is reduced by co- 
incident detections and follow-up observations, stacking 
analyses for nearby SNe could be useful (c.f. GRBs 31]). 

The detection of quasi-thermal neutrinos would prove 
the existence of relativistic neutron outflows, which can- 
not be probed by photon observations. We could get 
precious insights into magnetic acceleration mechanisms, 
the physics of PNSs, and clues to nucleosynthesis. Multi- 
wavelength studies at radio, optical and x-rays are also 
relevant to test our picture and constrain -B<jip and P. 
Detecting these neutrinos would also be useful in reveal- 
ing the SN-GRB connection. Indeed, rapidly-rotating 
and strongly-magnetized PNSs have been proposed as 
the central engine of successful jets leading to GRBs in 
compact stellar progenitors. A larger fraction of PNSs 
would have less extreme -Bdi p and P, leading to failed 
GRBs, hypernovae and perhaps super-luminous SNe. 

Importantly, quasi-thermal neutrinos, which we have 
discussed, do not require uncertain cosmic-ray ion accel- 
eration mechanisms, leading to non-thermal, high-energy 
neutrinos [29|, [3(j. For example, wakefield-like accelera- 
tion in the PNS wind could generate UHECRs and re- 
sulting PeV-EeV neutrinos via interactions with the stel- 
lar material and SN photons 29] . Note that mechanisms 
producing those neutrinos are completely different, which 
may occur only at t ;§> Uu n (where M » Mgj). 

The optical depth for the vN interaction is small. 
The near-spherical wind interacts with the shocked stel- 
lar material, and we have ~ 0.1(M sn /10M Q )i?- 2 10 at 
GeV, where R sn is the SN ejecta radius. The jet may 
interact with the pre-shocked envelope, and we have 
only - 10~ 4 at GeV if p ~ 1 g cm" 3 at 10 10 5 cm. 
However, studying neutrino oscillations in the McV- 
GeV range may be interesting. For example, the mat- 
ter effect is relevant since the resonance happens at « 
0.22 gem" 3 ( J B l/ /GeV)~ 1 (Am 2 /7.59xlO- 5 eV 2 ), and we 
can use it to probe the density profile or check neutrino 
properties given the progenitor structure (c.f. 32, 33j|). 
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